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Gravitropism, the slow reorientation of plant growth in response to gravity, is a major 
determinant of the form and posture of land plants. Recently a universal model of 
shoot gravitropism, the AC model, was presented, in which the dynamics of the tropic 
movement is only determined by the conflicting controls of (1) graviception that tends 
to curve the plants toward the vertical, and (2) proprioception that tends to keep the 
stem straight. This model was found to be valid for many species and over two orders 
of magnitude of organ size. However, the motor of the movement, the elongation, was 
purposely neglected in the AC model. If growth effects are to be taken into account, it is 
necessary to consider the material derivative, i.e., the rate of change of curvature bound 
to expanding and convected organ elements. Here we show that it is possible to rewrite 
the material equation of curvature in a compact simplified form that directly expresses the 
curvature variation as a function of the median elongation and of the distribution of the 
differential growth. By using this extended model, called the ACE model, growth is found 
to have two main destabilizing effects on the tropic movement: (1 ) passive orientation drift, 
which occurs when a curved element elongates without differential growth, and (2) fixed 
curvature, when an element leaves the elongation zone and is no longer able to actively 
change its curvature. By comparing the AC and ACE models to experiments, these two 
effects are found to be negligible. Our results show that the simplified AC mode can 
be used to analyze gravitropism and posture control in actively elongating plant organs 
without significant information loss. 
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INTRODUCTION 

Plant gravitropism is the capacity of plants to reorient themselves 
according to the gravitational field. The control over these move- 
ments involves several types of perception. In shoots, gravity is 
perceived by specialized cells all along the growth zone (Moulia 
and Fournier, 2009; Morita, 2010). This perception is then con- 
verted into a change in curvature of the organ powered by the 
differential growth mediated by enhanced lateral transport of the 
plant hormone auxin (Silk, 1984; Moulia and Fournier, 2009). 
Graviception has been an intensive field of study for more than a 
century (Sachs, 1888; Moulia and Fournier, 2009; Morita, 2010). 
However, it is only recently that kinematic studies have shown the 
universal importance of the straightening movement of organs 
where the sensing of curvature is local, e.g., aerial shoots (Bastien 
et al., 2013). This form of growth has been described in a sim- 
ple model, called the AC model. The AC model revealed that the 
control of the gravitropic movement involves not only gravicep- 
tion, but also proprioception of the organ curvature. For small 
curvatures, this model is written 



9C(s, t) 
dt 



-BsinA(s, r) - yC(s,t) 



where t is time, s is the curvilinear abscissa along the organ 
defined from the base s = 0 to the apex s = L(f), L(t) is the 
total length of the organ, A(s, f) is the angle of the material 
element at position s and time f with respect to the vertical, 
and C(s, f) = 3^ is the curvature of this material element 
(Figure 1). The proprioceptive term, —yC(s, t), tends to keep 
the plants straight. The graviceptive term — ft sinA(s, t) follows 
the sine law of gravitropism (Moulia and Fournier, 2009; Bastien 
et al., 2013) and tends to align the plant with the direction of 
the gravity field. When organ alignment approaches the verti- 
cal, near A(s, t) = 0, the graviceptive term can be linearized, 
ft sin A(s, t) ~ BA(s, t), yielding the most simple form of the AC 
model 

dC(s,t) 

= -pA(s, t) - yC(5, 0 (2) 

at 

The graviproprioceptive AC model (Equation 2) is both minimal 
and efficient in describing the observed universal pattern of cur- 
vature (Bastien et al., 2013). Equation (2) has been proposed as a 
way to phenotype gravitropism in shoot organs of plants, through 
the use of the dimensionless number B = ftL sz /y, where L gz is 
(1) the length of the zone that senses gravity and curvature and is 
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FIGURE 1 | Geometric description of a rod-like plant organ. (A) The arc 

length s is defined along the median line of the organ with s = 0 referring 
to the base and s = L the apex. The total length of the organ L can be 
broken down into the length of the growth zone L gz , which is moving with 
the apex, and Lf, the length of the zone which is not growing. This fixed 
portion elongates with time as elements leave the active growth zone. A(s) 
is the local orientation with respect to the vertical. (B) To simplify this, only 
two parts of the organ are considered to be elongating, the red part, at a 
position defined at its base s = Si, and the green part, at a position defined 
at its base s = S2- During elongation, the curvatures of these two parts are 
modified. Due to elongation, their distance relative to each other is 
modified. The position s= S2 no longer refers to the green element. To 
account for the variation in the material element, the derivative must follow 
the material element. 



able to undergo active bending. The number B quantifies the ratio 
between graviceptive and proprioceptive terms and entirely con- 
trols both the transient dynamics and the final steady-state shape 
of the gravitropic organ. When B was defined initially, the sub- 
script gz was coined for "growth zone," as the gravitropic response 
usually occurs within the growth zone (Morita, 2010). However, it 
is important to note that in the AC model the motor of the move- 
ment, i.e., the growth, was not explicitly specified. The movement 
was assumed to take place without any mean elongation so that 
the total length of the organ, L, and of the growth zone, L gz , 
remained constant. 

To move on from these assumptions to a more realistic 
description, it is important to extend the AC model by intro- 
ducing the motor of change of curvature, i.e., the differential 
growth across the stem, and the mean elongation (as plant cells 
in the primary growth zone do not shrink actively). To take 
into account the effects of growth, it is necessary to consider 
the change in curvature following the material points and not 
just the external shape even if this shape is stationary (Silk and 
Erickson, 1978; Chavarria-Krauser, 2006; Moulia and Fournier, 
2009; Merret et al, 2010) (Figure 1). In addition, through the 
elongation of the cells, the distance between material elements is 
continuously modified, and thus the orientation of the elements 
may change accordingly (Silk and Erickson, 1978). 

Here we investigate the consequences of including the mean 
elongation and differential elongation in modeling the control 
of the gravitropic movement of shoots. First, we extend the AC 
model to incorporate the effects of growth. As plants regulate 



curvature variation through differential growth, it would be use- 
ful to be able to express the visible output, the variation in 
curvature, as a function of differential growth. In order to prop- 
erly identify the effects of growth a simplified compact expression 
of this relationship is derived here. Second, the striking features 
of the AC model including elongation are explored. Two cases of 
elongation are studied, when the whole organ is elongating and 
when only the subapical part is growing. The destabilizing effects 
of these two types of growth are identified. The conditions for 
which the basic graviproprioceptive AC model (without growth) 
(Bastien et al., 2013) can account for the behavior of growing 
organs are discussed. Comparison with experimental data will 
give an estimation of the validity of these approximations. 

RESULTS 

KINEMATIC RELATIONSHIP BETWEEN CURVATURE VARIATION AND 
GROWTH 

When modeling an organ a given segment of the organ is approx- 
imated by a cylinder with a radius R (Figure 2). This is a 
sufficient approximation for the description of subapical plant 
movement like elongation and bending but not for apical growth. 
Gravitropic movement has been shown experimentally to lay 
within a single vertical plane defined by the initial orientation of 
the organ and the vertical. The curvature C of the organ thus lays 
within this plane (Darwin and Darwin, 1897; Correll and Kiss, 
2008; Moulia and Fournier, 2009). Within this plane both sides 
of the organ are growing with relative elongation growth rates, 
designated, respectively as ei and k 2 . From the simple geometri- 
cal model in Figure 2, the temporal evolution of the curvature 
of a material element with differential growth rate is given by 



DC(s, t) 
Dt 



1 

2R 



(1 



C(s, t) 2 R 2 )(e 2 



ei) 



(3) 



assuming there is no shear growth in plants [a plane cross section 
remains plane despite bending, fulfilling the Navier Bernouilli 
specification of the rod mechanics theory (Doghri, 2000), (see 
Supplementary Material for more details)]. As an organ expands, 
segments are moved along the organ relative to each other. 
So the material derivative comoving with each element of the 
organ has to be taken into account (see Silk, 1984; Moulia and 
Fournier, 2009 for more details). With respect to the curvilin- 
ear abscissa along the median line of the organ from the base 
to the apex (Figure 1), this material derivative is defined by 



DC(s, t) dC(s, t) 3C(s, f) 

= h v(s, f) 

Dt dt 3s 



(4) 



where v(s, f) is the velocity of the growth-induced displacement 
of each element along the organ. This is strictly equivalent to the 
expression given by Silk (1984), the first formula to express the 
relationship between curvature variation and differential growth 
(Supplementary Material). Another equation has been shown to 
be equivalent to the Silk equation up to the second order in 
C(s, t)R (Chavarria-Krauser, 2006). These three equations are 
thus equivalent, at least, for small curvatures. But the formulation 
of Equation (3) is advantageous here as it expresses the curvature 
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FIGURE 2 | A small element of an organ is considered as a cylinder of 
radius R with the upper side of length <Ssi and the lower side of length 

SS2- The sides are growing with a relative elongation rates of ei and <=2, 
respectively. After time f the curvature of the organ has changed. 



variation directly as the difference in elongation on each side of 
the organ. If i\ + 62 is not zero, i.e., if the organ undergoes net 
elongation, then we can separate the effects of the mean elon- 
gation rate of the segment, and of the relative distribution of 
differential growth across the organ as 



DC(s, f) 1 



Dt 



R 



(l - C(s, t) 2 R 2 )E(s, f)A(s, t) 



(5) 



DC(s, t)R 
Dt 



E(s, f)A(s, f) 



(8) 



Equation (8) expresses how the rate of change in curvature is pro- 
duced by the combination of the mean elongation rate and the 
relative amount of growth distributed to the differential growth. 
It thus expresses the motor role of elongation growth in the tropic 
movement. There is another element of growth that should be 
taken into account when considering tropic reactions. In many 
linear organs, like hypocotyls or primary roots, elongation is lim- 
ited to a zone just below the apex (Figure 1) (Silk, 1984; Gendreau 
et al., 1997; van der Weele et al, 2003). All the points at a distance 
greater than L gz have stopped growing, L(f) — s > L gz , where 1(f) 
is the total length of the organ from the apex at time t . As growth 
is the motor of curvature variation (see Equation 8), only the 
elements inside the elongation zone are able to undergo tropic 
bending. When a material element has been pushed out of the 
elongation zone and ceases elongation, E(s > 1(f) — L gz , f) = 0, 
its curvature cannot be modified any more (see Figure 1 in Moulia 
and Fournier, 2009). In most real cases E(s, f) displays a skewed 
bell curve, but as a simple starting point, the distribution of 
E(s, f) can be assumed to be a discontinuous step function that is 
stationary in time E(s < L(t) — L gz , f) = £0 (Figure 3). This sim- 
plification is expected to capture the main feature of the growth 
zone, that it first increases in size, L(t), and then has a con- 
stant length, L gz . The curvature variation is thus defined at every 
position s and at any time f as 



s < 1(f) - L. 



DC(s, t)R 
Dt 



(9) 



DC(s, t)R 

s > 1(f) - L gz ^— = £ 0 A(s, f) (10) 



where E(s, t) is the relative elemental elongation rate on the 
median line defined as the mean of the elongation rates on each 
side of the organ 



£(5, f) 



61 (s, f) + e 2 (s,f) 



(6) 



and A (5, f) is a dimensionless differential growth term which 
expresses the amount of growth distributed to the differential 
growth 



h(s,t)-€i(s,t) 

A(s, f) = 

61 (s, t) + e 2 (s, t) 



(7) 



It is important to note that E(s, f) and A (5, f) can be defined inde- 
pendently. If there is no shrinkage of the organ, €i(s, f) > 0 and 
£2(5, t) > 0, then I A I < 1. The characteristic length can be given 
by the radius of the organ R. By the same token, the characteris- 
tic time can be given by the inverse of the relative elongation rate 
of the median line E(s, f) , which gives the time for doubling 
the length [at constant £(5, f)]. As long as the radius of curvature 
C(s, t)~ l remains large compared to the radius R of the organ, i.e., 
CR 1, the quadratic prefactor in Equation (5) can be neglected, 
so the equation can be simplified 



So as the length of the organ L(f) increases over time, the length 
of the portion of the organ which can no longer curve, Lr{f), 
increases too. Another factor to consider is that there are dif- 
ferent phases of growth during organ development. In the first 
phase, when the organ is very young it elongates along its entire 
length. The length increases exponentially during this "exponen- 
tial growth phase." We can define the length of the growth zone 
as L gz . Some tissues then start to differentiate and cease extension. 
In this second phase, when (I > L gz ), only the subapical part of 
the organ is elongating, the length of growth zone L gz is almost 
constant, and the growth rate remains steady and maximal. The 
length of the whole organ increases approximately linearly with 
time during this "linear growth phase." In a third phase, which 
will not been studied here, the growth rate and the length of the 
growth zone decline until growth stops completely. 

CONTROL OF DIFFERENTIAL GROWTH BY GRAVICEPTI0N 

How does perception drive movement by modulating differential 
growth? When there is no differential growth, when A = 0, cur- 
vature can not be modified even if the organ is elongating, E > 0 
(see Equation 8). Furthermore, as there is no shrinkage of the 
organ, for a given A the relative elongation rate of the median 
can not modify the sign of the curvature variation, i.e., its direc- 
tion. Therefore we can assume that the perception affects A(s, f) 
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FIGURE 3 | The elongation rate of an organ. (A) When the organ is 


shorter than the length of the growth zone. 


L < L gz , it is expected that the 


whole organ is elongating. For simplicity, it is considered that the 


elongation rate is constant along the organ Eo- (B) If the length of the organ 


is longer than the elongation zone, L > L gz , 


the elongation rate is equal to 0 


outside the growth zone. 







and not E(s, t). Note that this is consistent with molecular studies 
of the distribution of the plant hormone auxin, which enhances 
cell elongation during the tropic response (Philippar et al., 1999). 
The total flux of auxin drives the mean elongation rate E while its 
transverse redistribution in tropisms drives the relative differen- 
tial growth (Silk, 1984). Keeping the graviproprioceptive equation 
of the AC model, but assuming that it only drives the distribution 
of differential growth A(s, t), yields 

A (5, t) = -MO, 0 - yC(s, t)R (ii) 

The tilde symbols on top of the parameters ft and y in 
Equation (11) is to indicate that the dimensions of these sensi- 
tivities have changed compared to the "non-growing" AC model 
(Equation 2) (Bastien et al., 2013) as the terms for differential 
growth in Equation (8) and 11 were broken down. The gravicep- 
tive sensitivity ft and the proprioceptive sensitivity y are now both 
dimensionless so that ft = RR/E and y = y/E. By doing this, it 
is now possible to directly compare the graviceptive and proprio- 
ceptive parameters of two different plants, independently of their 
size and mean elongation rate. The dimensionless number B can 
be expressed as: 

B=£f (12) 
yR 

Finally, from Equations (8, 11), the ACE model is defined by 

DC(s,t)R / ~ \ 

^=£(5,f)(-/iA(s,0-yC(U)fi) (13) 

Note that two plants with the same diameter R and the same sen- 
sitivities p and y can nevertheless behave differently at a given 
time t. For example, if the plants are elongating at different rates, 
the plant with the higher growth rate will converge to the vertical 
first. 



SOURCES OF EXPERIMENTAL DATA 

The experimental data used to assess the ACE model are those 
described in the study from Bastien et al. (2013). Experiments 
on etiolated wheat coleoptiles (Triticum aestivum cv. Recital) were 
conducted in growth cabinets. Experiments on bean hypocotyls 
(Phaseolus vulgaris), sunflower hypocotyls (Helianthus annuus), 
pea epicotyls (Pisum sativum), tomato stems (Solanum lycoper- 
sicum), chili stems (Capsicum annuum), raspberry canes (Rubus 
ideaus), carnation inflorescences (Dianthus caryophyllus) and 
Arabidopsis thaliana inflorescences (ecotype ColO) were con- 
ducted in controlled temperature greenhouses. Plants were grown 
until a chosen developmental stage of the organ of interest (e.g., 
until the beginning of inflorescence flowering for Arabidopsis). 
They were then tilted and clamped so the organ to be observed 
was horizontal so A(s = 0, t) = Ao = tt/2 for all f under con- 
stant environmental conditions in the dark (to prevent pho- 
totropism). The numbers of replicates were 30 for wheat, 15 for 
Arabidopsis, and 5 for all the other species. Data from similar 
experiments on Impatiens glandilufera stems by Pfeffer (Bastien 
et al., 2013) and on poplar trunks (Populus deltoides x nigra 
cv 14551) by Coutand et al. (2007) were also reprocessed and 
used. Simplified morphometric characteristics were measured 
following the methods described in Bastien et al. (2013). To take 
into account the two types of growth distribution (Figure 3), an 
effective length L e jf is defined as L gz when L > L gz and the ini- 
tial length of the organ Lq when L gz < L. The effective length 
L e jj was estimated by superimposing the first and last kine- 
matics images. The effective length of the organ is defined as 
the distance from the point where the plant starts to curve 
to the apex of the plant in the first image. The convergence 
length L c was obtained by plotting the local inclination angle 
A(s) along the organ beginning from the curved zone at the 
steady state. The angle A(s) was then fitted with the exponential 
A(s) = A 0 e- s/Lc . 

Here we first analyze the destabilizing effects of growth on the 
tropic movements and on the final shape of the organ, namely 
fixed curvature and passive orientation drift. We then analyze the 
behavior of the ACE model in the two cases of growth distri- 
bution along the organ, exponential growth when L < L gz and 
linear growth when L > L gz . Special attention will be paid to 
the situations in which growth effects may be neglected, when 
the simpler AC model may provide a sufficiently good esti- 
mate of the behavior of the more complete but complex ACE 
model. 

FIXED CURVATURE 

During the primary growth of plant organs, steady linear growth 
is usually achieved (in the absence of environmental stresses) in 
the first few days of organ development. Elongation then becomes 
restricted to a subapical zone whose length is L gz . Values of L gz and 
Eq remain almost constant for a few days (Figure 3B) (Silk, 1984). 
At this stage, the organ is only able to modify its own curvature 
from within the subapical growth zone (e.g., s > L(t) — L gz ). The 
length of fixed matured tissues convected out of the growth zone 
is given by 

L f (t)=L(t)-L gz (14) 
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If a material element that is convected out of the growth zone is 
curved, its curvature can no longer be modified and will remain a 
fixed part of the final shape. 

PASSIVE ORIENTATION DRIFT 

When there is longitudinal expansion (E(s, t) > 0) with no dif- 
ferential growth (A(s, f) = 0) (Figure 4), there is no variation 
in curvature j§ = 0. However, as the length of the element can 
increase due to elongation, the orientation of the two ends of the 
element might be modified. The difference in orientation A'(s, i) 
across a short segment Ss is given by A'(s, t) = C(s, t)8s. The vari- 
ation in orientation dA'(s, t)/dt of this segment over time can be 
expressed as (see Supplementary Material) 



dA'(s, t) 
It 



E 0 (A(s,f) + C(s, t)R) 



Ss 
R 



(15) 



This equation only applies to an independent element of infinites- 
imal length Ss, but not to the whole organ. As the perception and 
response are local (Bastien et al., 2013), the constraint on the reg- 
ulation of the shape must take place at the scale of each element 
independently of the other segments. 

When A(s, f) = 0 and the organ is curved, the orientation 
between the two ends of a material element is modified propor- 
tionally to its curvature. To achieve an active control over the 
movement, the differential growth term E(s, t)A(s, t) must be 
larger than the term for passive orientation drift, E(s, t)C(s, t)R, 
yielding 

|A(s,t)l > \C(s,t)R\ (16) 
| - /3A(s, t) - yC(s, t)R\ > |C(s, t)R\ (17) 

This process must be independent of the orientation so that the 
plant remains capable of active control even for small angles of 
perturbation, i.e., when A is close to 0 



y\C(s,t)R\ > \C(s,t)R\ 

y > 1 VA, C 



(18) 
(19) 



For an active tropic control to be possible, the proprioceptive term 
must overreach the effect of the passive orientation drift. The con- 
dition for this in an elongating system reduces simply to the value 
of the proprioceptive term y being greater than 1 . This means that 
the typical time for convergence to the steady state under the pro- 
prioceptive drive y~ l should be smaller than the typical time 
for expansion growth £ 0 . 



EXPONENTIAL GROWTH 

We can now consider the kind of growth that is observed in very 
young organs when the length of the organ increases exponen- 
tially over time (Figure 1). In this case the relative elongation rate 
E(s, t) is homogenous all along the organ and is homogenous in 
time (E(s, t) = Eg is considered to be constant). As seen previ- 
ously with passive orientation drift, according to Equation (19) 
the system cannot converge if y < 1 . Spatial oscillations of cur- 
vature along the organ are amplified over time. When the ori- 
entation of the organ tip reaches the vertical, the organ is still 
curved, and the effect of passive orientation drift spreads due to 
continued growth (see Figure 5A and Movie 1). Conversely when 
y > I the system can converge to a steady state (see Figure 5B 
and Movie 2). To maintain a steady shape despite the contin- 
uous disturbance from passive orientation drift due to growth, 
posture must be constantly regulated. This contrasts with the AC 
model for non-elongating organs where it is sufficient to balance 
graviception with proprioception. 

The steady state of the ACE model undergoing exponential 
growth is given by (see Supplementary Material for more details) 



A(s, f)=A 0 r(K) 



i v 



h- 



IPs 
R 



(20) 




8s 



FIGURE 4 | An element of a curved organ is elongating with the same 
growth rate i on each side. The increment of length does not affect the 
curvature of the organ but the final orientation of the element is changed. 



Bessel functions of the first kind /„ are oscillating functions, and 
r(x), the gamma function, is an extension of the factorial func- 
tion x\ to real numbers. As there are two Bessel functions, the 
steady state can stabilize some spatial undulations to produce a 
wavy shape (see Supplementary Material Figure 1). Such oscilla- 
tions can be avoided only if the first zero of the Bessel function 
is displaced toward the tip of the organ. One way to do this is 
to increase the order of the Bessel function, i.e., by increasing 
n = y — 1 in Equation (20) (so the first zero of ]„(x) is displaced 
to higher x values for increasing values of n). Only high y can 
avoid undulations in the steady state. As found previously for the 
control of passive orientation drift, the value of the proprioceptive 
sensitivity is central to avoiding undulations in the steady state 
shape, and hence to achieving proper gravitropism (Bastien et al., 
2013). 

CONDITIONS FOR NEGLECTING EXPONENTIAL GROWTH 

When does the simple AC model provide a good approximation 
for the behavior of the more complete but complex ACE model 
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FIGURE 5 | Simulations of gravitropic growth. The color (from blue to 
red) codes for the elapsed simulation time. (A) Solution of the ACE model 
during exponential growth when B = 10, y = 0.1 and Lq/R = 100. The 
simulated organ does not reach a steady state. The size of the curved zone 
is expanded by growth and the organ cannot regulate its posture (Movie 1). 
(B) Solution of the ACE model during exponential growth when 6=10, 
y = 10 and Lo/R = 100. The simulated organ reaches a steady state even if 
the organ is elongating (Movie 2). (C) Solution of the ACE model during 
subapical growth when B = 10, y — 0.1, and L gz /R = 100. The black curve 
is the part of the organ that is outside the growth zone where the curvature 
cannot be modified. As the simulated organ reaches the vertical, 
oscillations are fixed on the final shape due to the elements convected 
outside the growth zone (Movie 3). (D) Solution of the ACE model during 
subapical growth when B = 10, y = 10, and L gz /R = 100. The black curve 
is the part of the organ that is outside the growth where the curvature 
cannot be modified. The simulated organ reaches a steady state before the 
elements of the organ are convected outside of the growth zone. No 
oscillations are fixed on the final shape (Movie 4). 



for exponential growth? (See Supplementary Material for a math- 
ematical justification of the approximation of the steady states of 
the two models). 

We need to define "good approximation." As the AC and 
ACE models are intended to explain experimental results, a good 
approximation will be one where the difference between the AC 
and the ACE models is smaller than the experimental error in 
the "best" experiments. In our experiments the error in measure- 
ment of the orientation of the organ is of the order of 0.05 radians 
(Bastien et al., 2013). If the absolute difference between the steady 
state of the ACE model (Equation 20) and the steady state of the 
AC model (A(s) = A 0 e _Bs/i « z ), called /x, 



li= AqT(y) 



Bys 



i v 

2 



Jy-1 



Bys 



A 0 e~ Bs/I « z 



(21) 



is less than 0.05, it is not possible to distinguish between 
the two models. Numerically \i < 0.05 if y > 6.2 (Figure 2 



Supplementary Material). This approximate limit must, however, 
be expressed with measurable parameters from experiments if we 
are to test the validity of the approximation. As stated in the intro- 
duction, it is known that initiation of gravitropism in an organ 
results in an inhibition of growth on one side of the organ with a 
corresponding increase on the other side (Cosgrove, 1990; Myers 
et al, 1995; Tokiwa et al., 2006). According to Equation (7), this 
leads to |A| ~ 1. According to Equation (11), at the beginning 
of the tropic response, the distribution of differential growth is 
dominated by graviception so |A| ~ j8|Arj| ~ 1. However, this 
expression is only valid for small angles. Gravisensing at any 
angle was shown to follow a sine law (see Equation 1). Following 
the same argument as previously, it can be shown that |A| ~ 
B\ sin A| ~ 1. As | sin A| < 1, this yields R ~ 1. 

With /3 ~ 1 and Equation (12), a simple relationship between 
the proprioceptive term y, the length of the growth zone L gz and 
B can be stated as 



Y 



BR 



(22) 



The approximate limit defined for y can now be rewritten using 
measurable parameters B and L gz (Bastien et al, 2013) 



R 



> 6.2B 



(23) 



Equation (23) thus describes a condition for which the AC model 
gives an approximation of the steady state of the ACE model that 
lies within the range of experimental error. 

SUBAPICAL ELONGATION AND LINEAR GROWTH 

Growth is exponential only at the very start of organ develop- 
ment. Later the length of the growth zone is always finite, a generic 
feature of plant growth. Even in tiny coleoptiles, the length of the 
coleoptile soon overreaches L gz , and non-growing tissues accu- 
mulate at the base (Silk et al., 1989; Peters and Tomos, 2000; 
Walter et al, 2002). As the elements are convected out of the 
growth zone, the plant loses the ability to modify its curvature. 
Some plants overshoot the vertical during the gravitropic move- 
ment and display transient spatial oscillations (Bastien et al., 
2013). If elements are convected outside of the growth zone while 
overshooting the vertical, spatial oscillations could remain fixed 
on the final shape. Due to the complexity of the ACE model in 
this case, the steady state shape cannot be worked out analytically. 
However, the behavior of the model can be studied by numer- 
ical simulations. Figures 5C,D and Movies 3 and 4 show how 
the final steady-state shape (including fixed curvature) may vary 
depending on whether or not transient curvature becomes fixed. 
As all of the plants studied experimentally so far exhibited smooth 
convergence to the vertical without fixed oscillations of curva- 
ture (Bastien et al, 2013), it is important to clarify the conditions 
under which curvature undulation does not become fixed in the 
ACE model as a way of checking the validity of the modeling. 

CONDITIONS TO AVOID THE FIXATION OF CURVATURE 

If the organ has time to converge to a steady state before a ver- 
tical element is convected out of the growth zone, no spatial 
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oscillations are fixed on the final shape as there will be no inflec- 
tion point within the zone of fixed curvature. To assess whether 
spatial oscillations are fixed into the final shape, the time for a 
vertical element to be converted out of the growth zone, TV, has 
to be compared with the time to converge to the steady state, T c . 
If T c JS> TV, then a vertical element is convected outside of the 
growth zone while the organ is still regulating its posture. The 
organ cannot avoid fixation of transient over- curvatures during 
the overshoots of the vertical. However, if T c TV, when a ver- 
tical element is convected outside of the growth zone, the organ 
has already converged to its steady state, and this non-undulating 
steady-state shape is then fixed. 

The convergence time can be approximated from the proprio- 
ceptive term that dominates when approaching convergence 



1 



(24) 



To express TV, it is necessary to know the length of the elements of 
the organ that have left the elongation zone and lost their ability 
to curve. After time f the length of the organ is 



L(t) = L gz (l+Et) 



(25) 



By Equation (14), the length of fixed "mature" tissues convected 
out of the growth zone is given by 



L f (t) = 1(f) 



L gz Eot 



(26) 



To estimate the time required to fix a vertical element, we can con- 
sider an extreme case. If the whole plant is curving at the same rate 
all along the organ, depending only on the graviceptive term at the 
initial angle An, it will have a maximal bending rate. Assuming 
that this maximal rate is kept constant until the tip first reaches 
the vertical, a lower limit of the real time to reach the vertical 
can be estimated and should give the most stringent condition 
for neglecting fixed curvature. Thus the curvature of an element 
that leaves the growth zone at time r is given by 



C(L f (t), t)R = -PA Q E Q t 



(27) 



Integration over the fixed elements gives the angle reached by an 
element outside the growth zone 



A(L f (t), t)=A 0 - PAoE 0 t 



L f (t) 



(28) 



The first vertical element outside the growth zone is given simply 
by A(Lf(Tf), TV) = 0. From Equations (28, 26), the minimal time 
for a vertical element to leave the growth zone is given by 



2R 
fiL gz 



(29) 



The ratio between the time to converge to the steady state and the 



time to fix a vertical element, b£ 



^ £ , should be as low as 



possible for straightening to start and a straight vertical shape to 
be established. This defines a condition for proper straightening, 

IJET 



y V £ 5F' anc ^ tne higher the value the better the straightening. 
Equation (22) expresses this approximate limit as a function of 
parameters that can be measured experimentally, B and L gz where 



B < 



2L 



R 



(30) 



The dynamics of the organ is then confined close to the dynamics 
of the non-elongating model (with negligible curvature fixed in 
the mature basal zone). However, sub-apical growth restrains the 
range of values for the parameter y that make such straightening 
behavior possible. 

EXPERIMENTAL VALIDATION 

The value of B was found experimentally to range from 0 to 10, 
with most of the observed values ranging from 2 to 5 (Figure 6). 
Results in Figure 6 show that most of the plants tested fulfilled 
the conditions in which passive orientation drift can be neglected 
(Equation 23). Therefore in the range of experimental values of 
B, it is not possible to distinguish between the ACE and the AC 
models. It follows that the solution of the AC model without elon- 
gation is a good approximation of the ACE model. However, the 
growth of some wheat coleoptiles did not fulfil these conditions 
(Equation 23). It has already been noted that these coleoptiles 
oscillated less during the gravitropic movement than predicted 
by the AC model (Bastien et al., 2013). This behavior cannot 
be accounted for using the ACE model either, as the effects of 
growth identified in this paper are destabilizing effects that do not 
reduce oscillation. It is then likely that the effects of growth are 
negligible in the wheat coleoptile and that second-order mecha- 
nisms of regulation are involved. Overall the AC model that does 
not explicitly include elongation (Bastien et al., 2013) remains 
an excellent first-order approximation for an elongating organ. 
This is despite the effects of passive orientation drift and the fact 
that the observed final shapes in real plants are in fart dynam- 
ical steady states in which the organ is continuously expanding 
(described more precisely in the ACE model). 

The effects of fixed curvature can be neglected if B is in good 
agreement with the approximate bounds given by Equation (30). 
This condition can be experimentally tested by comparing L gz 
with B (Figure 6). All the plants of the different species and in 
all experiments were found to match the approximate lower limit 
given by relationship 30. Even if the length of the organ L is greater 
than the length of the growth zone L gz , the AC model is still a good 
approximation of the ACE model. 

DISCUSSION-CONCLUSION 

The ACE model developed here is an attempt to combine two 
branches of knowledge about gravitropism. Is it possible to 
describe the role of differential growth mediated by enhanced 
lateral auxin transport (Silk, 1984; Moulia and Fournier, 2009; 
Morita, 2010) in the context of the combined control of gravi- 
ception and proprioception formalized in the AC model (Bastien 
et al, 2013)? 
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FIGURE 6 | B = L eff /L c as a function of the effective length of the 
organ L eff /R. Each point denotes an individual plant. The red line 
represents B = ^/L^/R. All the points are under this line B < y/L B ff/R, 
showing that the destabilizing effects of growth on the whole gravitropic 
movement are negligible. The green line accounts for B = L Bff /6.2R. Most 
of the points are under this line, showing that it is not possible to 
distinguish the steady state of the ACE model from the AC model. 



The ACE model was built by extending the AC model with 
a simplified kinematic equation expressing how growth is the 
motor for the curvature variation in gravitropism (Equation 8). 
This equation directly relates the rate of curvature variation 
attached to a material element (an infinitesimal segment of the 
organ) ^ to the mean elongation E and the relative differen- 
tial growth rate A between each side of the organ. As in the AC 
model, the perception process combining graviception and pro- 
prioception drives the movement, but it is now assumed that this 
motor acts by controling the differential growth rate A (while the 
mean elongation rate remains constant). Introducing the motor 
of curvature variation in the model (Equation 8) has helped to 
identify two main effects of growth on the gravitropic response of 
shoots, and their consequences on achieving proper control over 
the tropic movement. 

First, the modification of the curvature can only occur when 
the relative differential growth is not equal to zero A / 0. Despite 
appearing trivial, this result contradicts the classical view of 
growth. It might be intuitively accepted that elongation of a 
curved cylindrical organ, without differential growth, decreases 
the curvature of this organ. Statements of this kind of intuition 
can be traced back at least as far as (Jost and Gibson, 1907), 
but the kinematics have not been demonstrated before. Here we 
show that if a material element is elongating with no differential 
growth, the curvature of this element is conserved. But the size of 
this element does increase. Therefore the orientation of its distal 
end, resulting from the integral of the curvature over the length of 
the element, is modified, due to curvature spreading over elonga- 
tion growth. This effect has been called "passive orientation drift." 
The consequence of such passive orientation drift is that the ori- 
entation of the growing curved organ must be regulated all the 
time for tropic control to be effective. Elongating organs can thus 



only achieve a purely dynamical tropic steady-state shape through 
a dynamical balance between the passive orientation drift and the 
active differential growth driven by graviproprioceptive control. 

Even when the shape remains stationary, the material ele- 
ments are "flowing" and their shapes must be actively regulated 
at all times. This is a new example of the maintenance of "steady 
shape from changing cells" (Silk and Erickson, 1978) which has 
been primarily described for the maintenance of the hook of 
some hypocotyls during their etiolation. More precisely, we have 
shown that to overcome passive orientation drift and reach a 
tropic steady-state shape, an approximate limit can be defined 
that relates the length of the growth zone L gz and its radius R 
to the bending number B where B represents the driving of the 
motion through the combined action of graviception and propri- 
oception all along the growth zone. Thus only organs with proper 
sets of values for the growth zone aspect ratio, L gz /R, and for 
the bending number B can control posture and convergent tropic 
movement. For example, organs with a high value for the bending 
number B will only reach a tropic steady state if they also have a 
long growth zone, L gz . 

The second effect of growth on the gravitropic response of 
shoots has been called the "fixed curvature." Again, this starts 
from an apparently trivial observation. All the elements outside 
the growth zone are no longer able to actively change their cur- 
vature or orientation [at least until secondary cambial growth 
begins (see Moulia and Fournier, 2009)]. But the consequence is 
not trivial and the ACE model is needed to elucidate this aspect 
of tropic growth. We showed that the organ must converge to the 
steady state shape fast enough so that only a negligible amount of 
transient curvatures remains fixed at the base. Experiments have 
shown that plants mostly exhibit a monotonic decrease in ori- 
entation from the tilted base to the almost vertical apex (Bastien 
et al., 2013). This is a good indicator that no curvature has become 
fixed. To avoid the fixing of transient curvatures, a new approxi- 
mate limit has been defined that relates the length of the growth 
zone, L gz , and the bending number B and thus the ratio between 
the graviceptive and the proprioceptive sensitivities. 

By revisiting the experimental data on gravitropic movements 
from several very different species and over two order of magni- 
tude in organ size, we have shown that all these plant movements 
are within the approximate bounds we defined by the model- 
ing. In other words, all the species and organ types studied so 
far have avoided fixed curvature and overcome passive orienta- 
tion drift (Bastien et al., 2013). Although gravitropism is thought 
to be adaptive (Moulia et al., 2006), there is very little evidence 
to support this. Our results may be the first evidence of converg- 
ing selection for phenotypic sets of (L gz , R, B) values that bring 
about efficient posture control. More direct studies of population 
genetics would be a way of investigating this hypothesis. 

The approximation of the ACE model with the AC model 
does not mean that growth can always be neglected. The tim- 
ing of the tropic movement is still fixed primarily by the mean 
relative rate of elongation growth. There are other aspects that 
could be included in modeling, especially the mechanical aspects 
powering the motor of differential growth (Moulia and Fournier, 
2009) or conditioning the mechanical stability of the straight- 
ened and erected shape. The tropic movement of the plant organs 
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with secondary growth in girth such as tree trunks and branches 
remains also to be studied. The motor for the curvature change in 
this case is not differential growth but differential shrinkages or 
expansion in the wood (Almeras and Fournier, 2009; Moulia and 
Fournier, 2009). 

The fact that plants were shown to fulfil the conditions 
required to overcome the two destabilizing effects of growth on 
plant movement has an overriding practical aspect. In these con- 
ditions, the original AC model is a good approximation for the 
more general but complex ACE model. This means that gravit- 
ropism and posture control can be characterized phenotypically 
by the measuring the dimensionless parameter B based only on 
the AC model. 
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